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Abstract 

The effects of random sources and sinks on the reaction kinet- 
ics have been investigated in the framework of field theoretic model, 
which is the result of Doi "second quantization" of corresponding mas- 
ter equation . It has been demonstrated that they do not affect strong 
impact of density fluctuations on the single-species annihilation reac- 
tion in two dimension. However random source nad sinks generate 
some new terms which have to be added into equation describing Gri- 
bov process. 

1 Introduction 

Reaction kinetics - main subject of our analysis belongs to the spreading pro- 
cesses, which attract much attention in last time due to their frequent occur- 
rence in physical, chemical, biological, ecological and sociological systems. 
Generally, such systems are characterised by large fluctuations of physical 
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quatities therefore it must be described by the methods of statistical physics. 
At the same time appropriate quantities, which give full information about 
system, must be introduced, e.g. probability distribution functions. For sys- 
tems with enormously large number of particles (e.g. gas, plasma etc.) they 
satisfy known chain of evolution equations. Reaction kinetics, which can be 
afilled to the system with large number of particles, operates with processes, 
where the number of particles is not conserved and probability distibution 
functions for fluctuating random quantities satisfy the master equations. Due 
to elegant and robust technique [1] these equations can be rewritten in the 
terminology of a "quantum theory" and after that by means of technique of 
functional integration can be formulated as the effective field-theoretic mod- 
els with a given action and defined random fields. These models describe 
mumerous processes with the same universality classes. It suggests situa- 
tion with various field- theoretic models of critical dynamics obtained from 
Langevin equation[2]. The Langevin equation is a stochastic equation with 
random forcing, which enters to the equation as additive or multiplicative fac- 
tor, compensates dissipative losses and sinks and, hence, ensures steady state 
of dynamics under consideration. In reaction kinetics the random sources and 
sinks, in fact, reflect real physical situation, in which during chemical reac- 
tion of follow up species, they can appear or disapear due to uncontroled 
random interaction with bath. , e.g. due to active chemical radicals. Theo- 
retical description of such a sources and sinks is not easy task. Their direct 
inclusion as defined additive regular terms and noises into Langevin equa- 
tion is artificial and their physical interpretation is dubious. There exists 
another physically more acceptable waynamely, introduction of new terms, 
responsible for interactions with bath, into corresponding master equation. 
In present paper we propose specific realization of general approach elebo- 
rated in [I] and we add the random sources and sinks to the master equation 
for annihilation equation. At the beginning we introduce the simplest terms 
responsible for the sources and sinks. Unfortunately, this simplest variant 
leads to non-conservation of particles and we have no possibility to com- 
pare it with standard Langevin approach. To correct this situation we add 
more complex sources, which keep the number of species, and technically 
can be compared with the effects produced by noices presented in Langevin 
equation. 

This paper is organized as follows. In section [2] we formulate general 
formalism for master equation for annihilation reaction with random sources 
and sinks. We recall basic features of Doi "second quantization" formalism 
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and we construct basic field theoretic dynamical action fuctional. Section [3] 
is devoted to the analysis of possible situations, which can occur, when we 
consider different space dimensions and canonical dimensions of master fields. 
Section @] is devoted to discussion of the results and concluding remarks. 

2 Random sources and sinks in the master 
equation 

We will consider the annihilation reaction A + A — > in a random drift field 
in a more general setup than in [3] . To this end we introduce random sources 
and sinks for the particles, e.g., to maintain a steady state in the system. 
In most cases this is carried out by including an additive noise term in the 
Langevin equation of the stochastic process. Since our analysis is based on 
the master equation, this is not quite appropriate here. Unfortunately, there 
is no unique way to introduce random sources in the master equation corre- 
sponding to the random noise of the mean-field (Langevin) description. We 
use the simplest choice, described in detail, e.g., in jl], which is tantamount 
to amending the problem by reactions A — > X and Y — > A, where X and 
Y stand for particle baths of the sink and the source, respectively. In a 
homogeneous system these reactions give rise to the master equations 



where P(t, n) is the probability to find n particles at the time instant t in the 
system. The ellipsis in ([[]) stands for the terms describing the annihilation 
reaction, diffusion and advection in the system. In (JTJ) fi + and /i_ are the 
reaction constants of the creation and annihilation reactions, respectively. 
The transition rate has been chosen proportional to the particle number 
n, which seems the quite natural and also preserves the empty state as an 
absorbing state. In the transition rate for creation process V is the volume 
of the (for the time being) homogeneous system and will be important in 
passing to the continuum limit of the inhomogeneous system. We recall that 
the master equation ([!]) gives rise to the reaction-rate equation 




= n+V[P{t,n-l)-P(t,n)] 

+ n-[(n + l)P{t,n+l)-nP(t,n)]... 



dt 



(1) 
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where (n) is the mean particle number. 

The idea of the Doi approach [T] is to rewrite the set of master equations 
for probability distributions of a stochastic problem in the form of a single 
kinetic equation for a state vector incorporating all probabilistic information 
about the system constructed in a suitable Fock space. The kinetic equation 
is defined by the Liouville operator acting in the Fock space and generated 
by the set of master equations. Although the basic procedure has been 
thoroughly exposed in the literature, the introduction of random sources and 
sinks of particles in the master equation has specific features which should be 
presented in detail. Therefore, let us briefly recall the basic quantities and 
relations of the Doi approach. For simplicity, consider probabilities P(t, n) 
to find n particles at the time instant t on a fixed lattice site. The spatial 
dependence may then be described by labeling the particle number by the 
coordinates of the lattice and introducing necessary sums and products over 
the lattice sites. 

Basis vectors of the Fock space are defined by the usual annihilation and 
creation operators a, d + and the vacuum vector | 0): 

a|0) = 0, a + \n) = |n + l) , (3) 
[ a, d + ] = aa + — a + a = I , (4) 

with the normalization (n\m) = n\S nm . 

The state vector containing all information about the system is defined 

as 

oo 

|$) = ^P(t,n)|n). (5) 

n=0 

The set of master equations for a birth-death process may be cast in the form 
of a single evolution equation for the state vector (jSJ) without any explicit 
dependence on the occupation number: 

dl 

-^ = L(a + ,a)|$>. (6) 

Master equations (JT|) give rise to the following terms in the Liouville operator 

L g (a + , a) = n + V (a + — /)+//_(/ — d + ) a , (7) 

where I is the identity operator. The expectation value of any function F(n) 
of the random particle number 

oo 

<F(t)) = £>(n)P(f,n), (8) 

n=0 
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may be expressed in the form of the functional integral over the functions of 
time a and a: 

(F(t)) = J VaDaF N [(a(t) + l)a(t)]e Sl , (9) 

where F^(aa) is the normal form [5] of the operator F(a + a) and Si is the 
dynamic action 

oo 

St{a, a) = Jdt \-a(t)d t a{t) + fi + Va(t) - n_a(t)a(t)] . . . (10) 
o 

Only the generic time- derivative term and terms brought about by the ran- 
dom source model are expressed here explicitly, while the ellipsis stands for 
terms corresponding to other reactions and initial conditions. 

Let the transition rates /z± be random functions uncorrelated in time with 
a probability distribution given in terms of the moments (//") = E± >n . At 
this point we also generalize the treatment to the case of a spatially inho- 
mogeneous system and introduce a lattice subscript as the spatial argument, 
e.g. a(t) — > cii(t). In this case the volume V becomes the volume element 
attached to the lattice site. For simplicity, we replace the time integral with 
the integral sum f£°dt — > Y2 a ^ an d assume that the transition rates at 
each time instant and lattice site /i±,a,i are independent random variables. 
Then the average of the expectation value ([9]) over the distribution of random 
sources reduces to the calculation of the expectation value 

J_ J_( e ^ + • a ' iV ^°'' i ^' t ~ IJ '~ .o^a^n.iAi^ (11) 

a.i 

For each particular time instant and lattice this yields (we assume that the 
moments of fi± are the same for all a and i and omit labels for brevity) this 
gives rise to the usual cumulant expansion 

(e^ A *) = l + bAtE 1 + lE 2 (bAt) 2 + lE 3 (bAt) 3 + --- 

2 o 

_ e bAtE 1 +i(E 2 -E2)(bAt) 2 +±(E 3 -3EiE2+Ef)(bAt) 3 +- ^ 2 ) 

Here, b stands for either Vd or —da. Thus, for instance the average over /i + 
assumes the form 

J^^ e M+, a>i Va aj< A^ = ^^[AtE+rVa^ + ^E+i-El^iVaa^At) 2 } 
a,i 

x e E a E i [|(E+3-3E +1 E +2 +E| 1 )(y5a,iAt)3+..| ( 13 ) 
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In the continuum limit the function a Qji is replaced by the field ip + (t,x), 
whereas in the limit V — > the expression a a ^/V gives rise to the field 
ijj(t, x). The sum over a together with At gives rise to the time integral and 
the sum over i together with the volume element gives rise to the spatial 
integral ^ V — > J dx. In the first term of the exponential in (fT3l) we thus 
obtain 

Yl AtE +i v ^ -> E +1 dt rfx^ + (t,x). 

a i J J 

In the cumulants of second and higher order the continuum limit is not so 
obvious. We assume the simplest nontrivial distribution for fi±, in which 
only the variance term has a finite limit, when At — > and V — > 0, whereas 
the contributions of higher-order cumulants vanish, for instance 

(E +2 - E 2 +1 ) VAt a + , At , V ->• , (14) 
(£ +3 - 3E +1 E +2 + E 3 + ,) (VAt) 2 , At -> , V -> . (15) 

Therefore, the contribution of the average over fi + to the effective dynamic 
action assumes the form 

S + = Jdtjdx j£ +1 ^ + (t,x) + l -a + [^ + (t,x)] 2 | . (16) 

For the average over //_ a similar argument yields 

S_ = jdt jd* |-E_ 1 ^ + (t,x)^(t,x) + ia_ [^ + (t,x)^(t,x)] 2 | . (17) 

These contributions to the effective dynamic action may, of course, be gen- 
erated by suitably chosen normal distributions of fi±. 

This way of introduction of random sources and sinks has the annoying 
feature that it does not conserve the number of particles in the system. For 
a comparison with the treatment of this problem in the Langevin approach 
the random sources and sinks should be introduced in such a way that the 
particle number is conserved. 

The simplest way to effect this is to add to the random source a term 
proportional to the particle number, i.e. use the "reaction constant" fi + V + 
/ii + n instead of fi + V in the master equation. The source terms on the right- 
hand side of the master equation ([1]) in this case the assume the form 

= » + V[P(t,n-l)-P(t,n)) 

+ jj l+ [(n-l)P(t,n-l)-nP(t,n)}... (18) 
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The new part of the master equation corresponds to a branching process [3]. 

The added term gives rise to the following contribution to the Liouville 
operator 

Lg2{d> + , a) = (a + — /) a + a . (19) 

Performing the steps described above we arrive at the contribution to the 
dynamic action in the form 

s 1+ = Jdt Jd* (^ + + 1) v + \°i+^ +2 + 1) V} • (20) 

It is easy to see now that if we exclude the plain source (i.e. put E +1 = a + = 
0) and choose = E_i, then the empty state remains absorbing and the 
"mass term" oc disappears in the dynamic action and we arrive at the 
dynamic action of random sources and sinks 

S gc = Jdtjdx | J E 1+1 V +2 ^ + \o- (ip + if>) 2 

+ Ia 1+ ^ + > + + l)Vj (21) 

which conserves the average number of particles. 

The effects of the high-order terms are drastically different in the two 
cases amenable for a scaling analysis with the aid of the renormalization 
group. The time derivative term in the dynamic action 

S=- Jdt Jdxip + (t,x)d t i)(t,x) + ... 

must be dimensionless in order to have nontrivial dynamics. Therefore the to- 
tal scaling dimension of the number-density operator ip + {t, x)?/>(t, x) is equal 
to the dimension of space and thus positive. 

First, if the scaling dimension of the field ip + is equal to zero, d^+ = 0, 
then the dimension of the field ip is positive (more precisely, d^ = d) and 
the operator monomials in the second and third terms in ( I2T1) have the same 
scaling dimension. Since they are carrying the factor ip 2 , the their scaling 
dimension is larger than that of ip +2 ip. Therefore, the second and third 
terms in ff2T]) are IR irrelevant and should be discarded in the asymptotic 
analysis. Second, if the scaling dimensions of both fields are positive, then 
in the operator monomials in the second and third terms in ( 121]) there is 
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at least one "excessive" field factor in comparison with the first term which 
renders them irrelevant. Thus, in these cases the IR relevant dynamic action 
of random sources and sinks reduces to the single term 

S' gc = J dt J d^E 1+l i) +2 ^p , d f + = V d$+ > 0, d,p > . (22) 

Third, if the scaling dimension of the field '0 = 0, the scaling dimension of 
the field ip + is positive, terms with "excessive" powers of ip + are IR irrelevant 
and the starting point for the subsequent RG analysis is the source and sink 
action in the form 

2 



S'; c = dt dx{ + 7: (<t- + <ti+) h <h = • ( 23 ) 



3 Annihilation reaction A + A — > with ran- 
dom sources and sinks 

Let us now analyze the dynamic action of the diffusion-limited annihilation 
reaction A + A — > 



DC 







S 1 = - dt d*{ ij + d t ij - J D ^ + VV 



+ X D [2^+ + (^+) 2 ] ^ 2 }+n J dx^+(x,0) . (24) 

from the point of view of scaling behaviour sketched in section [3 

In the first case with d^+ = the nonlinear terms in action fl24l) are of 
equal scaling dimension. However, the source-sink part (1221) is linear in the 
field ip with positive scaling dimension in contrast to the quadratic in if) terms 
of (|24p . Therefore, the IR relevant interaction term in this case is ( I22p and 
the corresponding dynamic action is 

Sim = -J dtjd* {0 + d t 0> - ,Do^ + VV - 

+ n o y"rfx^ + (x,0). (25) 

However, this dynamic action does not bring about any graphs with closed 
loops of the density propagator and thus no density fluctuation effects on the 
asymptotic behaviour are anticipated. 
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In the second case with > and > the fourth-order term in ac- 
tion becomes irrelevant. Either of the remaining third-order terms alone 
does not generate loops, therefore density fluctuation effects are brought 
about only, when both fields have the same scaling dimension d^+ — — 
d/2. In this case the IR relevant dynamic action is 



oo 







Sua = ~ dt dx< ij + d t ij - D ij + V 2 4> 



+ 2X D Q ip + ip 2 - E 1+1 ij +2 ij^+n J rfx^+(x, 0) . (26) 

This is the dynamic action of the Gribov process [5], also known as the 
Reggeon model, subject to random advection. Effects of the random drift 
with the use of the Obukhov-Kraichnan compressible velocity field have been 
analyzed in [7]. 

In the third case with d^ = the fourth-order term in action (|24l) becomes 
irrelevant as well due to the positive dimension of the field By the same 
token, however, both terms of the source-sink action (|23|) are also irrelevant 
and we arrive at the IR relevant dynamic action 



oc 



Sim = -J dtj dx - A)^ + V^ + 2A A)^V 

+ n Jdxip + (x,0). (27) 

This dynamic action does not give rise to any density-fluctuation loops and 
thus does not predict any decay anomaly due to them. 

In summary, if the sources and sinks are chosen such that they conserve 
the mean number of particles in the system, then the anomalous scaling 
behaviour in the system is that of the Gribov process, if any. 

A different situation arises, if the plain source term is included in the 
analysis. Then there is the possibility that the system does not tend to the 
absorbing empty state but to an active state with a finite concentration of 
particles. In this case the starting point is the dynamic action with all the 
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terms quoted above, i.e. 

S = jf dtjdx^-^ + d t ^ + D ^ + V 2 ^-\ D [2^+ + (^+) 2 ] ^j 2 + E +1 ij + 

+ \*+ ty + f + E l+l ^ + (V + + 1) i> + \o^ +2 (V + + l)V " £-1^ 
+ ia_(^) 2 |+n o yrfx^ + (x,0). (28) 

The stationarity equation brought about by this dynamic action for the field 
■ip is (the stationary value ip + = as usual) 

dt^fj - A)V 2 ^ = -2A A^ 2 + E +1 + E 1+1 ^j - . (29) 

However, the action expanded around the stationary value brought about by 
this equation is rather complicated. To keep expressions simple, continue to 
consider the case E i+i = E_ x . Then the re-expanded action is 

S = J^dtJ d^-i) + d t ip + Do^+V 2 ^ - V8y/E +1 X D 7p + 7p 

V2X D 4A A) 4A A) 2 1 



2A 0j D 4A A) AqDo 



^E +1 X D <T 1+ 1p + 1p „ / (7i+ (7_\ +2/2 

+ a 1+ + ai+ ^ 2 } + n y dx ^ + (x, 0) . (30) 

In the critical limit E'+i — >■ 0. Since it is the expectation value of a nonnega- 
tive random quantity //+, the variance cr + vanishes as well. In the vicinity of 
the critical point we keep only the leading E +x and a + putting them equal 
to zero in terms, where they are subleading. This simplifies the action a bit 
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S = J dtj d*{-i) + d t i) + D ip + V 2 ip - V8y/E +1 X D ^ + ilj 

/ WggT o+\ 2 E +1 a 1+ ^ E +1 a 1+ ^ + 2 

V V2X^ 2X D 4X D lV v 

+ V2^E +1 Xj\a 1+ jj +3 ^ | y/E +1 X D Q a 1+ tjj +4 ilj 2 \ D ib + ib 2 

Ao-Do \/2 XqD 

( \ n i ai + i a ~\ l+ 2 I 2 , 1+3/2 , cr 1+ lp +4 ip 2 \ 
+ [-XqDq + — + —J ip + ip +a 1+ ifj + if) + 1 

+ n dx^ + (x,0) . (31) 



Dimensional analysis of the canonical dimensions then yields the following 
cases. In the nonlinear parts without the critical parameters E+\ and a + the 
previous arguments hold, but in terms having powers of these parameters 
as coefficients the positive scaling dimensions of them must be taken into 
account. The free-field part of the action (j3~Ti) suggests that the canonical 
dimension of E + i is four. The canonical dimension of cr+, in fact, remains a 
free parameter. 

Proceeding in the same manner as above, we arrive at the following effec- 
tive actions for the IR scaling limit. In the first case with d^+ = the third 
and fourth powers of if> + and independent of if> or first order in if) due to the 
coefficients proportional to E+\ or its square root are irrelevant compared 
with terms oc ip +2 in action (l3"Tj) . Nonlinear in if> terms are irrelevant against 
the linear terms due to positive dimension of if). Therefore, the IR effective 
action in this case is 

S = J dt J dx {-ip + d t ip - ip + V[vV>] + D ip + V 2 ifj -2V2 y/E +1 X D if> + if) 

+ (^^i + ^^ + ^^l+no/rfx^O). (32) 

Again, the only remaining interaction terms does not bring about loops, 
although here we have a nontrivial correlation function of the field if). 

In the second case with d^+ > and > higher powers than the 
leading corrections to the free- field action of both fields are irrelevant. This 
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argument leaves us with the dynamic action 



S = J dtjd*{ -i) + d t i) - tp + V[vf] + D tfj + V 2 tp -2^2 y/E +1 A A) ^ 
+ n cix^ + (x,0) . (33) 



Contrary to the case discussed above, here the interaction term —2 Ao-Do ip + ip 2 
generates loops alone due to the presence of the correlation function of the 
field ip. Therefore, two effective actions with nontrivial fluctuation contribu- 
tions are possible. 

a) d^p+ > d^p. To keep the correlation function of the field i[) for the loops, 
the variance cr+ must have a dimension less than that of y/E + \. This yields 
the effective action 

S = dt j d-x {--0 + SfcV - ^ + VW + LW + V 2 V - 2 V2 V E +1 X D ^ + tlj 

+ ^^ +2 -2Ao J Do^V} + ^y'rfx^ + (x,0) (34) 

with the critical dimension depending on the scaling dimension of cr + in the 
spirit of the description of tricritical scaling behaviour [2]. The model is 
logarithmic at six dimensions, however, because apart from the coefficient of 
the oc ip +2 the action is that of critical dynamics of the <p 3 model. 

b) d^+ = dj/j = d/2. Both third-order terms are relevant and the effective 
action is basically fl33|) . In this case the dimension of a + is larger than that 
of y/E + i and for simplicity we omit cr+. Thus, the effective dynamic action 
may be written as 

S = dt Jdx{-^ + d t ip -^ + V[v^] + D ijj + V 2 ij -2V2y/E +1 \ D iP + ij 
t ~ l x/7r ^ +2 + £-i^ +2 ^-2A A)^ + ^ 2 i + ™o /dx^ + (x,0). (35) 



Note that this is a dynamic action describing the Gribov process with a 
random source independent of the active agent density. That the rate of 
change of the density due to the random sink is proportional to a power of 
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density is a natural assumption. The assumption that the rate of change of 
the density due to the random source is proportional to a power of density is 
not natural. Therefore, the dynamic action fl35|) possibly predicts a critical 
behaviour of the Gribov process different from that discussed in the literature. 

In the third case with > and = we arrive at the effective action 
AMD- 

The analysis of scaling dimensions shows that we may actually lift most 
of the restrictions on the probability distribution of the transition rates of the 
type ([H]) and (fT5|) . Indeed, even if the higher order cumulants are finite, the 
scaling dimensions of corresponding terms in the dynamic action grow with 
the order of the cumulant with the exception of the case, when the transition 
rate is independent of the agent density. 

4 Conclusion 

We have investigated possible effects of random sources and sinks on annihi- 
lation reaction. Contrary to frequently used approach, in which the sources 
and sinks are introduced directly to the Langevin equation, we have added 
it directly to the master equation, where their physical sense is clear. On 
the basis of dimensional analysis we have constructed effective actions, which 
are a starting point for RG analysis of critical behaviour of systems under 
consideration. Most of the dynamic effective actions obtained are trivial 
from the point of view of critical behaviour of the model, because infrared 
relevant interaction terms do not generate loops diagrams. We can conclude 
that for the study of influence of density fluctuations on annihilation reac- 
tion A + A — > 0, which are relevant only for dimension two or less, random 
sources and sinks are unimportant. Rest non-trivial cases lead to the action 
describing Gribov process. However, in this manner we obtain a source term 
in dynamic action without agent density. This is cardinally different from 
that generated by Langevin equation, therefore RG analysis can lead to the 
new effects in critical behaviour, when sources and sinks asymptotically van- 
ish. We need analyse dependence of scaling functions on sinks and sources 
parameters in infrared limit. It suggests the situation, which takes place in 
phase transitions, where statistical correlations of order parameter depends 
on "mass" (deviation of temperature from Currie one) and dependence of 
their scaling functions on '"mass" have to be investigated. This analysis we 
intend to carry out in near future. 



13 



The work was supported by VEGA grant 0173 of Slovak Academy of Sciences, 
and by Centre of Excellency for Nanofluid of IEP SAS. This article was also 
created by implementation of the Cooperative phenomena and phase transi- 
tions in nanosystems with perspective utilization in nano- and biotechnology 
projects No 26220120021 and No 26220120033. Funding for the operational 
research and development program was provided by the European Regional 
Development Fund. 

References 

[1] M. Doi, J. Phys. A: Math. Gen. 9 (1976) 1465, 1479 . 

[2] A. N. Vasil'ev, The Field Theoretic Renormalization Group in Critical 
Behavior Theory and Stochastic Dynamics (Chapman&Hall/ CRS Press 
Company New York, 2004). 

[3] M. Hnatich, J.Honkonen, Phys. Rev. E 61, 3904 (2000). 

[4] N. G. van Kampen, Stochastic processes in physics and chemistry 
(North-Holland, Amsterdam, 1984). 

[5] A.N. Vasiliev, Functional Methods in Quantum Field Theory and Sta- 
tistical Physics (Gordon and Breach, Amsterdam, 1998). 

[6] V.N. Gribov, Zh. Eksp. Teor. Fiz. 53, 654 (1967). 

[7] N. V. Antonov, A. S. Kapustin, J. Phys. A: Math. Theor. 43, 405001 
(2010). 



14 



